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A well-established method to deal with highly correlated systems is based on the expansion 
of the Green's function in terms of spectral moments. In the context of the Composite Operator 
Method one approximation is proposed: a set of n composite fields is assumed as fundamental basis 
and the dynamics is considered up to the order s. The resulting Green's function has a sn-pole 
structure. The truncation of the hierarchy of the equations of motion is made at the s-th order and 
the first s — 1 equations are treated exactly. A theorem, which rules the conservation of the spectral 
moments, is presented. The procedure is applied to the Hubbard model and a recurrence relation 
for the calculation of its electronic spectral moments is derived. 

(N 

£4 , I. INTRODUCTION 

The discovery of new materials with unusual properties^ requires the formulation of new schemes of calculation 
capable to catch the essential physics of correlations among particles. One powerful tool is the Green's function 
formalism. Most of the properties of a many-body system can be fully described in terms of the thermal single-particle 
' Green's function. For finite systems the spectral function p(k, lu) can always be written as an infinite weighted sum 
, of delta functions and the retarded Green's functiona can be expressed as 

o 
o 



g% u )= = V » (i.i) 



/<Q (k,q/)~ trW(k) 

^ lu - Ei(k) + ir] 

where £7j(k) and <jW(k) are the energy spectrum and the spectral density function, respectively, of the i-th elementary 
excitation. In the way of developing approximate methods, one line of thinking is to approximate the Green's function 
by a sum on a finite number of excitations (pole approximation) . The main question that arises is how the poles must 
be chosen (i.e., which excitations do describe the physics we are interested to analyze?^. Let us consider a Hamiltonian 
H = H[tp], where tp denotes a field operator. In the spectral density approach (SDA)u one truncates the sum in (1.1) 
by picking up the first s terms (s-pole approximation). The unknown functions -Ei(k) and erW(k) are fixed by the set 
"jS ■ of 2s nonlinear equations 



m (k) i 



>(k)=$>*(k)a»(k) (1.2) 
where m^ fc ^(k), the moments of the spectral density, are calculated by means of the equations of motion: 

O 

m« (k) = F.T. \{id/dt x ) k < [ V {x), (y)] ± >] (1.3) 

• I L J t x —ty 

5^ ■ Eq. (1.2) is derived by making use of the exact relation (conservation of spectral moments) 



W(k) = / dLup(k,u)iu k (1.4) 



The symbol [• • -]± denotes the anticommutator (+) or the commutator (— ) according to the statistics the field ip 
obeys; the symbol < • ■ ■ > denotes the thermal average on the grand canonical ensemble; F.T. indicates the Fourier 
transform. This procedure, which is an expansion of the Green's function in terms of the first 2s spectral moments, 
gives a good description in the high-energy region but certainly fails in the low-energy region, where one needs s to 
be very large. 
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When we consider systems where the correlation is very high the procedure above illustrated may not be adequate. 
Due to the strong interactions the original particles are not observable; new excitations appear and it is more convenient 
to describe the system in terms of new field operators. Collective behaviors in forms of bound states, resonances 
and diffused modes emerge as the physical fields. Although some of them are not stable excitations, they give 
considerable contributions to the physical processes and therefore it is sometimes convenient to promote them to the 
role of well-defined quasi-particle excitations. The choice of new fundamental particles, whose properties have to be 
self-consistently determined by dynamics, symmetries and boundary conditions, becomes relevant. The necessity of 
developing a formulation to treat composite excitations has seen remarkable developments in the many-body theory- iii 
the form of an assortment of techniques that may be termed composite particle methods. The projection methodofl, 
the method of equation of motiorjj, the slave boson methodQ, the spectral density approach^, the composite operator 
method (COMpu are all along this line. In the COM language let us write H = Hj + Hjj, where Hj is chosen 
according to the physics of the problem under investigation, and let ip = {^pi} [i = 1, • • • , n] be a complete set of 
eigenoperators of Hj 

i^(x) = Mx), Hj] = e/HV)V(*) (1.5) 

The fields ip n are generally constructed as functionals of the original field tp and are called composite fields. Then, Hn 
will be taken into account by some approximate methods that will use ip as operatorial basis. Various approximations 
have been used to estimate the contribution coming from Hn. One method, the n— pole approximation, linearizes 
the full equation of motion 

d - 
i-g-Mx) = [ip(x), H] = e(-W)tp(x) + Sj(x) w e(-iS/)tp(x) (1.6) 

where the eigenvalue or energy matrix e is self-consistently calculated by means of the equation 

e(k)M(°)(k) = AfW(k) (1.7) 
Af( fc )(k) is the generalized spectral moment, given by the n x n matrix 

M«(k) = F.T. \{id/dt x ) k < [iP(x),^(y)} ± >] _ (1.8) 

The decomposition in (1.6) is chosen in such a way that the higher-order source dj(x) does not have a projection 
on the fundamental basis {tpi} ■ In this approximation the retarded "single-particle" Green's function G(k, ui) = 
F.T. < R[^{x)^{y)} > has a n-pole structure, where the number of poles is equal to the number of composite fields. 
It is clear that the properties of the Green's function are controlled by the dynamics and by the choice of the new 
basis. When this approximation is used the COM has some similarity with the SDA: both approaches correspond 
to expanding the Green's function in terms of spectral moments. However, there is a fundamental difference. In the 
SDA there is no freedom in choosing the spectral moments, which are completely fixed by the equations of motion. 
At the contrary, in the COM the spectral moments are mainly determined by the choice of the operatorial basis. The 
idea is that in the framework of approximating the Green's function by a finite number of excitations, a proper choice 
of the basis might give a better convergence and a more accurate description of the low-energy region. Let us write 

M W (k) = (k) + <5A/W (k) (1.9) 

where 

M^(k) =e fc (k)M (0) (k) (1.10) 



fe-i 

,5Af( fc )(k) = ]T F.T.e h ~ m - x (id/dt x ) m < [Sj(x),^(y)} ± > t ^ ty (1.11) 

m— 1 

Af( fe )(k) corresponds to the spectral moment calculated in the n-pole approximation, where the higher-order field 
Sj(x) in the Heisenberg equation (1.6) is neglected. In Ref. [lC] we proved the following theorems. 

Theorem I Given a set of composite fields {ipi, i = 1, n}, satisfying linearized Heisenberg equations i(d/dt)i/; = 
eip, the approximate spectral moments M^^k), are conserved at any order. 
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Theorem II Given a set of composite fields {tpi, I = 1, ....,n} , if the subset {ipi, I = 1, ....,n — 1} satisfies linear 

2 + 1 

Heisenberg equations i-^ipi(x) = lip{~~ ^) q Pp( x ) i then the first 2(n — I + 1) spectral moments for the field ipi [that 

P =i 

is Mjp (k) ] are conserved. 

The Theorem II shows that the SDA corresponds to COM when a particular choice of the basic fields is considered 
and a pole approximation is used. In COM we have more flexibility, in the sense that there is no restriction on the 
choice of the basic fields. 

A different approximation has been proposed in Refs. |Tl| and [fll, where in the context of the t-J and Hubbard models 
the incoherent part of the Green's function has been taken into account by a two-site approximation in combined 
use of the resolvent method. Although this approximation gave excellent results, the numerical calculations are very 
heavy in the small temperature regime. A pole approximation is much more convenient for practical calculations and 
there are no restrictions on the range of physical parameters. As shown in a variety of works, the pole approximation 
gives a satisfactory account of the local and spectral features. Therefore, it is interesting to explore the possibility 
to improve the approximation of the COM by improving the n-pole approximation. In this article we propose a 
new approximation, the sn-pole approximation, where we combine two ingredients: the choice of the fields and the 
dynamics. The procedure is the following: (i) we fix n basic fields by some reasonable criterion based on the physical 
properties of the model; (ii) we take into account the dynamics up to the s-th order. In this scheme the Green's 
function has a sn-pole structure 

fz(v- Ei{k.) + ir\ 

where -Ej(k) and <jW(k) are expressed in terms of the 2s spectral moment matrices M^(k) [0 < k < 2s — 1], Let 
A( m '(k) [m = 1, . . . , s — 1] be matrices of rank n x n determined by the equations 

s-1 

M (s+k) (k) - A {m) (k)M {m+k \k) =0 [for 0<fc<s-I] (1.13) 

m=0 

Then, the main results of this work can be summarized in the following theorems. 

Theorem III Given a set of composite fields i = 1, . . . ,n}, satisfying linearized Heisenberg equations 

£) \ s s — 1 / rj \ ni 

l dl) ^)-E AW HV) ^ ^)=0 (1.14) 

' m=0 ^ ' 

the first 2s generalized spectral moments are conserved and the approximate spectral moments M^(k) are conserved 
at any order. 



Theorem IV Given a set of composite fields {ipi, I = 1, . . . ,n}, if the subset {ipi, I = 1, ,,..,n — 1} satisfies the 

L U 



Heisenberg equations (1.14), the first 2s(n — I + 1) spectral moments for the field ipi(x) [that is M,^(k)] are con- 



served. 

When a particular choice of the basic fields is considered, we will show that the sn-pole approximation is equivalent 
to SDA with s — > sn. 

II. EQUATIONS OF MOTION AND GREEN'S FUNCTION 

The composite field ip{x) satisfies the equation of motion 

i^{x) = [i>{x),H]=jV{x) (2.1) 

The source j^(x) is expressed, in general, in terms of higher-order composite fields with a dynamics determined by 
the Heisenberg equation 
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i^-/ 1 Hx) = \jW(x),H]=jW(x) (2.2) 

The process can be continued and, in general, we have for the 5-th derivative the equation 

Pi \ s 

i-J xb(x)=j^(x) [for s>l] (2.3) 

Let us rewrite Eq. (2.3) under the form 

Pi \ s / Pi \ m 

f at WW-E AM HV) 1>{x) = 8jM(x) (2.4) 

' m=0 ^ ' 

where AS m \— iV) are matrices of rank n determined by the equations 

<[5j^( x ,t),(-i^-) ^(y,t)]±>=0 [for < k < s - 1] (2.5) 



By substituting (2.4) into (2.5) we find that the matrices ^ m )(k) [0 < m < s — 1 ] are determined by the following 
system of equations 

s-l 

M( s+fe )(k) - ^ (m) ( k )^ (m+fe) ( k ) = [for 0<k<s-l] (2.6) 

The evaluation of the s matrices A^ m ^(k) requires the knowledge of the first 2s generalized spectral moments M^(k), 
that can be calculated by means of (1.8). By using the equation of motion (2.4) the retarded Green's function 
G(x,y) =< R[ip(x)tp^ (y)] > satisfies the equation 

E^(-^) (i-^-) G(x,y) 

k=0 \ xJ 

= E»* (ar) V* "*„) E A< m \-iV x )ll m - k \x,y) (2.7) 

fc=l ^ X ' m=k 

- <R[Sj^(x)^(y)} > 

where 

I (k \x,y)=<[j (k \x),^(y)} ± > t!c=tv (2.8) 

with the convention that 

A^(-iV) = -1 (2.9) 

Eq. (2.7) is an exact equation. We shall define our approximation, sn-pole approximation, by neglecting the term 
< R[5j( s \x)^(y)] > , with 5j^(x) fixed by the requirement of vanishing projections on the basis \j^ k \x); < k < 
s — 1] with the convection j^(x) = ip(x). Then, the solution of (2.7) in momentum space is 

G(k, lo) = —± £ ° (S,k) M"*" 1 ( 2 - 10 ) 

lo s - A( fe )(k)o; fc fe=i 

k=0 

where we defined 

s 

C (s > k) (k) = - Y A^(k)M (m - k \k) [for l<fc<s] (2.11) 

m—k 

Eq. (2.10) shows that G(k, uj) is expressed in terms of the first 2s generalized spectral moments. This expansion is 
dictated by the choice of the composite field and by the dynamics up to the s-th order. Expression (2.10) can be 
rewritten in the spectral form 
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G(k,u,) = £. 



<r«(k) 



^ w - £j(k) + ir? 
where the energy spectra Ei(k) are the roots of the equation 



g u/"a m (k) = 



(2.12) 



(2.13) 



m=0 



with coefficients a m (k) defined by 



s-l 

E 

fe=0 



g^DetK- E A«(k)u, fc } 



a m (k) = - 
a sn (k) = 1 

The spectral density functions a^> (k) are expressed as 



[for < m < sn — I] 



(2.14) 



sn — 1 



where we put 



CT(<) ( k ) = WfkTS^(k)AM (k ) 



&w(k)= n [^(k)-^(k)] 

3=^,3^ 



(2.15) 



(2.16) 



and the A^ m ^ (k) are n x n matrices, determined by the recurrence relation 



A«(k) = E A(f)(k)A( fc+s -P)(k) 



p=0 

s-l 



-EE a fe+s -p(k)A(" l )(k)M(" l -P- 1 )(k) 

p—0 m— p+1 



[for < k < sn - 1] 



(2.17) 



with the convention 



AM(k) = for m<£[0,sn-l] 
a m (k) = for m ^ [0, sn] 



(2.18) 



When k > sn — 2s the matrices (k) can be simply expressed in terms of the spectral moments by means of the 
relation 



i-k-l 



A( fe )(k)= «-- m (k)M(— fc - x - m )(k) 



[fc 



2s<k< 



1] 



(2.19) 



m=0 



In particular, we see that for n < 2 the relation (2.19) is valid for all values of s and for k varying in the all interval 
< k < sn- 1. 

Summarizing, the calculation of the Green's function in the sn-approximation requires the following steps: 

• (i) given a Hamiltonian, we choose a set {ip(x)} of n composite field; 

• (ii) we calculate the first 2s generalized spectral moments M^(k) by means of the formula (1.8); 

• (iii) we calculate the matrices j4/ m )(k) [0 < m < s — 1] by solving the linear system (2.6); 

• (iv) we calculate the "characteristic" coefficients a m (k) [0 < m < sn] by means of (2.14); 

• (v) the energy spectra Ei(k) [1 < i < sn] are calculated as the roots of the equation (2.13); 

• (vi) the matrices A^"^(k) [0 < m < sn — 1] are calculated by means of the recurrence relation (2.17); 

• (vi) the spectral functions a^(k) [1 < k < sn] are calculated by means of the expression (2.15). 

In closing this section, it is worth noting that the truncation the hierarchy of the equations of motion at the order 
s implies that the first s — 1 equations are treated exactly. 
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III. CONSERVATION OF THE SPECTRAL MOMENTS 



The spectral moments play an essential role in the formulation developed in the previous Section. It is important 

to verify to which extent the exact relation (1.4) is conserved in our approximation. In the sn-pole approximation 
the spectral function has the expression 

sn 

p(k, W )=2aW(k)5[ W -£7 i (k)] (3.1) 
»=i 

Let us denote by B^(k) the spectral moments calculated in the present approximation 

+?° 

B< fc >(k) = / 4jwV(k,w) = ^£f(k)ff (i) (k) (3.2) 

The point to discuss is the relation between B^(k) and the exact spectral moments M^(k) , calculated by means 
of (1.8). We shall study separately the two quantities. In Appendix A we prove that the matrices B^(k) satisfy the 
following relations 



(i)For n=l 



(k) = (k) for < k < 2s - 1 



s-l 



(ii) For n > 2 



SW(k) = - X) a m (k)£( fe - s+TO )(k) for k>s 

m=0 



BW(k) = M< fc )(k) for < fc < 2s - 1 

fi(fc) (k) = £ AW (k)B ( - k - s +P'> (k) for s < k < sn - 1 

sn — 1 

fi(fc)(k) = - ]T a m (k)£( fc - s "+ m )(k) forfc>sn 

m=0 



(3.3) 



(3.4) 



We see that the first 2s spectral moments are conserved. This is true for any element of the matrix M^(k), and 
therefore for all the fields tpi [1 < I < n]. In Appendix B we show that the spectral moments for fc > s can be written 
as 

M {k) (k) = M {k) (k) + <5M (fe) (k) (3.5) 
where iVfW(k) and 5M<- k \k) satisfy the recurrence relations 

s-l 

M«(k) = A ( - m \k)M < - m+k - s \k) [for fc>s] (3.6) 

s-l 

5M (k \k) = F< fc - S >(k) + Y A (m \k)SM^ n+k - s \k) [for fc > 2s] (3.7) 



m=0 



with 



F«(k) =F.T. < [(i^j Sj^(x),^(y)} ± > tx=ty (3.8) 

Use of time translational invariance and of the constrains (2.5) leads to 

<5M W (k) = [for < fc < 2s - 1] (3.9) 



G 



We note that (k) represents the spectral moment in the approximation considered, when the higher-order field 
Sj^ is neglected. Furthermore, in Appendix B wc show that the matrices M^(k) satisfy the recurrence relation 

sn— 1 

M {k \k) =-J2 a m {k)M {k ~ sn+m \k) [for k > sn] (3.10) 

By comparing (3.10) and (3.3)-(3.4) we see that for k > sn Af( fc )(k) and B^ k \k) satisfy the same recurrence relation. 
But, we have shown that for k < sn — 1 

B (k \k) = M {k \k) (3.11) 

therefore the two matrices coincide at any order. This shows an internal consistency of the approximation, in the sense 
that the exact relation (1.4) is preserved at any order when calculated on the same ground. The above discussion can 
be summarized in the Theorem III stated in Section I. 



IV. SPECIAL CHOICE OF THE BASIS 

Up to now the basic set ip(x) has not been specified. Let us consider the case where the n fields are chosen as 

(a \ s l+l 8-1 I / 8 \ m 

%) M*) = E^Hv-^O*) + EE< m) Hv) Ydi) Mx) (41) 
' p=l m=lp=l ^ ' 

where is a matrix of n rows and n + 1 columns and A 1 - 171 * 1 [1 < m < s — 1] are n x n matrices. That is, we choose 
the first field ipi (x) ; then, all the other fields tpi [2 < I < n] are determined by means of the equation 

/ r) \ s l — l / o \ m 

AfX^-wmx) = (i-j ^-iw-Ea-Uh^)^) 



-EEA ( -i P HV)( 

m-lp-l ^ 



ig- t ) M*) (4-2) 



It is clear that our basis satisfies the equation of motion 

(r) \ S S ~ * / fl \ m 

^)=E A(ro) HV)N^J Hx)+6j^(x) (4.3) 
' m=0 ^ ' 

where A^ is the square matrix obtained from A^ by taking the first n rows and n columns. The matrices A^ , 
< m < s — 1] and the vector Sj^ have the structure 

2(0) _ for 1 < I < n - 1 

A ^p " U /0r l+2<p<n+l 

A ( : } =0 for )f l ^ n - 1 < 4 " 4 ) 

l 'P J I + 1 < p < n 

Sj^ =0 for l<n-l 

We shall take the matrices A^ [1 < m < s — 1] in such a way to satisfy the constrain (2.8). We have seen in Section 
3 (cfr. Eq. (3.5)) that the spectral moments can be written as M^(k) — M^(k) + 5M^ (k) . For the special choice 
of the basis considered in this Section, it is possible to show that the elements of the matrices (k) , defined by 
Eq. (3.8), satisfy the following rules 



(k) = for I < n — 1 [for any value of p and k] 
Fi k ^(k) = fork<s(n-p+l)-l 



(4.5) 



Then, the elements of the matrix SM^(k) satisfy the following rule 

5m[i ] (k) = for k< 2s(n - I + 1) - 1 (4.6) 
By recalling the result B^ (k) = M( fe ) (k) , we have thus proved the Theorem IV, stated in Section I. 
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V. SOME EXAMPLES 



In order to illustrate the formulation let us consider some specific examples. 
Case n = 1 

The solution of the system (1.17) is 



A^(k)=J2(-^-) M^+'-^k) [for 0<k<s-l] (5.1) 

1=1 V-^W / k+l.l 

where D(k) is the s x s matrix 

/ M(°) ••• M( s -D \ 
I>(k)= : ... : (5.2) 
y M^" 1 ) • • • M( 2s " 2 ) y 

Then, the retarded Green's function is given by 

G M = ± fffi . (5.3) 

The energy spectra Ei(k) are determined as the roots of the equation (2.13), with characteristic coefficients a m (k) = 
_ ^4.( m )(k) . The spectral functions are given by (2.15) with 

s-fe-l 

A<*>(k) = - ^ A( s - m )(k)M( 5 - fe - 1 - m )(k) [.for < k < s - 1] (5.4) 

m=0 

This case is exactly equivalent to SDA and the previous formulas give a complete solution to the problem of solving 
the system (1.2) for any value of s. 
Case n = 2 

For n — 2 we introduce the basic field 

(55) 

and consider the particular case where the second field is taken according to the special choice (4.1). By noting that 
the equation of motion for the field ip2 {x) is 

(i£) k M*) = -[^(-^V)]- 1 £ A^(-iW)j[ k+m \x) (5.6) 

the generalized spectral moments are given by 

M^ik) =m(*>(k) 

M«(k) = M<*>(k) = -[^(k)]- 1 £ A^ ) (k) TO ( fe +")(k) 

m=0 V -'- 1 

Mi 2 fe) (k) = [Ag(k)]- 2 £ E^^ ) (k)4 P i ) (k)m( fc +™+rt(k) 

m— p— 

We see that in order to calculate the matrices A^ m>> [0 < m < s — 1], we need to know the spectral moments 
77j( fc )(k) , relative to the first field ipi(x), for < k < As — 1. Lengthy calculations show that the equation (2.6) can 
be expressed as 

2s 

^ a m (k)™( fc+m )(k) = [for < fc < 2s - 1] (5.8) 

m=0 

This result shows that the case (n — 2,s) is equivalent to the case (n = 1,2s), when the special choice (4.1) is 
considered. We can make the conjecture that for this particular choice the formulation depends only the product sn. 
That is, two approximations characterized by the pairs (m, Si) and (712, S2) are equivalent if n\Si = n 2 S2- 
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VI. THE HUBBARD MODEL 



As an application of the formulation let us consider the Hubbard model, described by the following Hamiltonian 

H = - l*8ii)J(i)cU) + U^ni(i) ni (i) (6.1) 

ij i 

The first term is the kinetic term which describes the motion of the electrons among the sites of the Bravais lattice 
described by the vectors Ri [i — (Ri,t)] ; c(i),c^(i) are annihilation and creation operators of electrons at site Ri 
in the spinor notation 

-(el) ct =(4 4) (6-2) 

and satisfy canonical anti-commutation relations: 

K(i),4,(j)} = S^Sij K(i),<v(j)} = {<£(*), c t,(j)} = (6.3) 

fi is the chemical potential, tij denotes the transfer integral and describes hopping between different sites. We fix the 
scale of the energy in such a way that t a = . In momentum space the hopping matrix can be written as 



't(k) (6.4) 



where iV is the number of the sites, d is the dimensionality of the system, ft is the volume of the unit cell in the 
coordinate space and f2s is the volume of the Brillouin zone. In addition to the band energy the model contains a 
term approximating the interaction among the electrons. In the simplest form of the Hubbard model the interaction 
is between electrons of opposite spin on the same lattice site; the strength of the interaction is described by the 
parameter U . n a (i) = c) a {i)c a {i) is the charge density of the electrons with spin a. The equation of motion for the 
electron field is 

i—c(i) = &(-iV)c(i) + Uri(i) (6.5) 

where 

6(-*V) = -/i + t(-iV) . . 

r,{i)=c^{i)c{i)c{i) ( 1 

In the framework of the sn-pole approximation we shall consider the case n = 2, s = 2 . Let us take as basis the 
two component field 



and let us consider two possible choices. 
I Choice 

= c(i) 



Mi) = -[^(-^V)]- 1 £ A^(-iV) (i£) m c(i) (6 - 8) 

m=0 

This field satisfies the equation of motion 

a \ 2 1 / o \ m 

^) V>«=£^ (m) Hv)(^j + (6-9) 

' m=0 ^ ' 

with the constrain 

l 

M (2+fc)( k ) _ A( m \\t)M( m+ V(k) = [0 < k < 1] (6.10) 

m=0 
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As shown in Section 5, by taking this choice the formulation (n=2, s=2) is equivalent to the case (n=l, s=4). Then, 
the retarded Green's function has the expression 

G ( k,c) = y ^ (6.11) 

The energy spectra are the roots of the equation 

4 

a m (k)cj m = (6.12) 

with coefficients determined by the equations 

4 

a m (k)m (fc+m) (k) = [0 < k < 3] (6.13) 

m=0 

The spectral functions are 

-j 3 m 

ffW(k) = WW) E ^ (m) (k)E s r m_1 (k) ap (k) (6.14) 

^ ' m=0 p=0 

We need to calculate the first 8 electronic spectral moments 

m^\k)=F.T. <{(i^r) c(i),ct(j)}> B . T . (6.15) 

The calculation of the electronic spectral moments can be done by means of a recurrence relation, as shown in 
Appendix D. 
II Choice 

This field satisfies the equation of motion 

! 1 / p. \ rn 

^(x)-Y J A^\-K)[i-\ <p(x)=6j(x) (6.17) 

' ' ' m=0 ^ ' 

with the constrain 

l 

M (2+fe) _ A {m) {k)M {m+k) (k) = [0 < k < 1] (6.18) 

m=0 

The spectral moments M^ k \k.) = F.T. < {(idij){i) / dti) k , V'(j)}ti=t :j > can be expressed in terms of the electronic 
moments (6.15) by means of the following relations 

M[ k) (k) = m( fe )(k) - 2d( fe )(k) + .g( fe )(k) 

M x ( 2 } (k) = M^\k) = rf( fe )(k) - ff (*)(k) (6.19) 
M 2 ( 2 fe) (k)=g( fe )(k) 

where 

(k) = U- 1 [m( fe+1 ) (k) - 6(k)m« (k)] . , 

.gW(k) = J7- 2 [m( fe+2 )(k) - 26(k)m( fe+1 )(k) +6 2 (k)m«(k)] ( J 

We need to calculate the first 6 electronic spectral moments. The retarded Green's function has a four-pole structure 
with energy spectra and spectral functions calculated as shown in Section 5. As stated in Section II, we note that 
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the approximation has been made on the equation of motion for second-order time derivative. No approximation has 

been done on the equations for the first -order time derivatives. In particular, the following equations 

if £(*) = 6(-iV)£(i) + t(-*V)»7(*) - 7r(i) , fi 2 , 

i|»?(i) = (C/- M )r?W+7T(i) 1 ' j 

where 

7r(i) = ^n^(i)[t(-iV)c(i)} + c(i)[t(-iV)J(i)]c(i) (6.22) 
are treated exactly in this approximation. 



VII. CONCLUSIONS 



Among various analytical methods developed in the last decades to deal with highly correlated electronic systems, 
the technique of expanding the Green's function in terms of a finite number of spectral moments has been rather 
successful in describing the properties of many correlated models. Along this line of thinking, the crucial point is the 
choice of the basis in terms of which the expansion is realized. A standard way is to choose the spectral moments 
of the original fields. However, in presence of strong correlation this choice may not provide the better basis for a 
perturbative approach and other basis may be considered. In particular, a set of composite field operators may be 
taken as the fundamental set in terms of which the Green's function formulation is realized. In this process there are 
two ingredients that should be considered: the choice of the basis and the dynamics. In this work, we have presented 
a scheme of calculations where both aspects are taken into account. Firstly, a set of n composite basic fields is taken; 
the choice of this set is made on the basis of physical criteria: the role played by the interaction and by the value 
of the external parameters, the phase we are interested to describe, and so on. Secondly, the dynamics is considered 
up to the s-th order: we consider time derivatives of the fundamental fields of increasing order and the truncation 
procedure is applied to the equation of motion for the time derivative of order s, with the advantage that all equations 
of motion up to order s — 1 are treated exactly. Two theorems have been proved about the conservation of the spectral 
moments in the approximation considered. Finally, we have applied the formulation to the case of the Hubbard models 
where some specific realizations have been considered. In Appendix D we have presented a recurrence relation for the 
calculation of the spectral moments relative to the Hubbard model. 
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APPENDIX A: CONSERVATION OF THE FIRST 2S SPECTRAL MOMENTS 

The spectral moments calculated in the sn -pole approximation have the expression 

sn 

B( fe )(k)=^^(k)aW(k) (Al) 
»=i 

By using the expression (2.15) for the spectral functions, we have 

sn— 1 

5 (fe) (k) = \ {m) (k)p( k+rn \k) (A2) 

where we have defined 

P W(v) = y!$M (A3) 

2 — 1 v 7 



11 



By observing that the energy spectra satisfy the equation 

sn 

a ro (k)Pf(k) = [for i = l,---sn] (A4) 

m=0 

the quantities p( fe ' (k) satisfy the recurrence relation 

sn 

^a m (k)p( fe+m )(k) = [for k>0] (A5) 

m— 

By simple algebraic relations we can show that 

pW(k) = (? {° r 0< fc < 5 -- 2 ( A6) 
y ' [1 /or fc = sn — 1 v y 

Then, Eq. (A2) can be written as 

sn — 1 

B (fc) (k) = A( m) (k)p( fe+ro )(k) [/or 0</c<sn-l] (A7) 

m— sn — 1 — 

When < k < 2s — 1 we have always sn — 1 — k > sn — 2s ; then, we can use Eq. (2.19) to obtain 

sn — 1 sn — m— 1 

B(fc)(k)= £ £ M( s "-" l - 1 -P)(k)a s „_p(k)p( fe + m )(k) 

m—s7i—l — k p— 

sn 

= £ M("0(k) £ a p (k)P(f+ fe " 1 -'")(k) 

m— p=sn— k+m \^ L ^ > ) 

k—1 sn 

= MW(k)a s „(k)p( s "- 1 )(k) + £ M( m )(k) £ a p (k)p(P+' £ - 1 - m )(k) 

m— p— sn — k-\-7n 

= M( fe )(k) [for < fc < 2s - 1] 

where Eq. (A6) has been used. Then, we have the result that in the present approximation the first 2s spectral 
moments are conserved. In closing this appendix we note that from the above discussion the following relations follow 
( i) For n = 1 



(ii) For n>2 



pW (k) = ikfW (k) for < k < 2s - 1 

pW(k) = - 5 a m (k)p(' £ - s + m )(k) /or fc > s 

m=0 



BW(k) = MW(k) for0<fc<2s-l 

s-1 

P( fe )(k) = £ A(P)(k)p(' £ - s +P)(k) for s < fc < sn - 1 

sn— 1 

P( fe )(k) =- £ a m (k)pC £ - s "+ m )(k) for k > sn 

m=0 



APPENDIX B: 

From the field equations (2.4) we have for k > s 

k s — 1 / ^ \ m+k — s / ^ \ k — s 



(A9) 



(A10) 



\ K S — ± / r\ \ m-\-K — s / p\ \ k — s 

V>(z) = £^ (m) HV)(^J ^) + (^J « w (») (Bi) 



m=0 

Then, for k > s the exact spectral moments satisfy the equation 
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s-l 



AfW(k) = A^(k)M^ m+k - s ^(k)+F^ k -^(k) 



m=0 



where we defined 



F (fc) (k) = F.T. < 



>t x =t v 



Use of time translational invariance and of the constrains (2.5) leads to 

F (k-s) ( k ) = o [for s <k<2s-l] 

Then, the spectral moments for k > s can be written as 

M (fe) (k) = M (fe) (k) +(5M (fe) (k) 
where M< fe )(k) and SM^(k) satisfy the recurrence relations 

s-l 

Af( fe )(k) = A^ n \k)M^ m+k - s \k) [for k>s] 

m=0 

s-l 

<5M W (k) = (k) + A(m) {k)5M( m+k - s ) (k) [for k > 2s] 



m=0 



with 



8M {k) (k) = [for < fc < 2s - 1] 



(B2) 

(B3) 

(B4) 
(B5) 

(B6) 
(B7) 
(B8) 



We note that M^{k) represents the spectral moment in the approximation considered, when the higher-order field 
Sj^ is neglected. In Appendix C we prove the following: 

Theorem V Given s matrices [0 < m < s — 1] of rank n, let us define the characteristic polynomial as 

s— 1 sn 

£>(w) = Det[-uj s + A^uj m ] = Y u m a m (B9) 

m— m— 

and let us consider the n x n matrices X^ m ^ [m > 0] determined by the following recurrence relation 

s-l 

X^ = J2 A^ n '>X ( - m+k -^ [for k>s] (BIO) 



m=0 



Then, the following relation holds 



X^ k+m ^a m = [for k > 0] 



(Bll) 



m=0 



for any choice of the matrices X^ [0 < m < s — 1]. 

By means of this theorem the matrices M^(k) satisfy the recurrence relation 

sn — 1 

M«(k) = -J2 a m {k)M { - k - sn+m \k) [for k > sn] 

m=0 



(B12) 
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APPENDIX C: PROOF OF THEOREM V 

Given s matrices [0 < m < s — 1] of rank n, we consider the matrix 



P(u) =lu s I-J2 A (m) Lo m (CI) 



m=0 

where w is a parameter, which can generally be complex. / is the unit matrix in n dimensions. Let us suppose that 
there exists the inverse matrix 

s-1 

Q(lu) = [lu s I - J2 A^u" 1 }- 1 (C2) 
This matrix can be always written as a rational function with respect to the parameter u) 

«"> - w> (C3 > 

The numerator B(ui) is a matrix whose elements are polynomials in u of degree s(n — 1) . The denominator ip(u) is 
the characteristic polynomial of degree sn 

s— 1 sn 

<p(u>) = Det[-uj s + A {m) iu m ] = ^ m «™ ( C4 ) 
Since Q{oj) is the inverse matrix of P(w) we have 

s-1 



[lu s I - J2 A {m ^ m ]Q{uo) = I (C5) 

m=0 

Use of Eq. (C3) leads to 

s-1 

[u s I - A (m) uj m }B(uj) = <p(w)I (C6) 

m=0 

The two members of this identity are polynomials of degree sn with matricial coefficients. In order that this identity 
be satisfied, the coefficients for each power uj m must be equal. In this identity we can replace the variable u m with a 
matrix X< m+fe ) [k > 0] 

Lu m -» X {m+k) (C7) 
If the matrices _X"( m ) are chosen in such to satisfy the recurrence relation 



X ( S +k) = j2 A (™) X ( m +k) ^ or k > j ( C8 ^ 

m=0 

then the left hand in (C6) vanishes, and so the second: 

sn 

- a mX {k+m) = [/or k > 0] (C9) 

m=0 

for any choice of the matrices X^ [0 < m < s — 1] . 
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APPENDIX D: SPECTRAL MOMENTS FOR THE HUBBARD MODEL 



We recall the definition of the electronic spectral moments 



,(fe) 



Qi)=F.T.<{jW(i),c(j)}> E . T . 



with 



(Dl) 



(D2) 



In order to calculate the anticommutators, we shall proceed in the following way. From the Heisenberg equation (6.5) 
we obtain for k > 1 



d 



fc-i 



(D3) 



with the convention that = c(i) . By recalling the definition of the Hubbard operator rj(i) [cfr. Eq. (6.6)] we 

have 



4) = E(m)E(T) (-^jf(H m - l) ^ (k - m) (i) 



m =o v 7 2=0 
The field equation (D3) can be written for k > 1 as 



(D4) 



= &HV)/ fe " 1 )(i) 
fe_1 / k — 1 

m=0 



rn 



E 



(-l) , ^ )t (i)^ m - ,) (i)i (fc - 1 - m) (i) 



(D5) 



We have a recurrence relation which allows us to calculate the source j^ k \i). Let us define the quantities 



E.T. 



9 k af3 (iJ) = {j { a k)i (i),4(j)}E.T. 



(D6) 



By means of the previous results, these quantities satisfy the following recurrence relations for k > 1 

f k a ^3) = b{-iV i )f k a0 \i,j) 



k-l 



+^E 



m=0 



k - 1 

m 



(D7) 



m=0 ^ ' 1=0 ^ 7 



with 



s^(i.i) = o (D9) 

Then, the spectral moments are given by 

m( fc )(k) = fe(k)m( fe - 1 )(k) + L/L( fe - 1 )(k) [for k > 1] (D10) 
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where 



4V 1) ( k ) = '£ V/ ' 1 

m=0 



with 



E ( "i ) (-l) l F.T. 



<j^(U m - l Hi)f^ 1 - m (i,j) > E . T . 
-<jf{i)Q-%j)j [ t 1 - m) {i) >E.T. 

+ <9 i P ^j)j ( r l) (^ 1 ' m) ^) >e.t. 



(Dll) 



»<°>(k) 



(D12) 
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